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ABSTRACT 


The Schwarz-Neumann alternating method is employed to obtain stress 
intensity solutions to two crack problems of practical importance: (1) a 

semi- inf inite elastic plate containing an edge crack which is subjected 
to concentrated normal and tangential forces, and (2) an elastic half- 
space containing a semicircular surface crack which is subjected to 
uniform opening pressure. The solution to the semicircular surface 
crack is seen to be a significant improvement over existing approximate 
solutions. Application of the alternating method to other crack prob- 
lems of current interest is briefly discussed. 
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ALTERNATING METHOD APPLIED TO EDGE AND 
SURFACE CRACK PROBLEMS 
by 

R. J. Hartranft and G. C. Sih 
Lehigh University 


1 • INTRODUCTION 

Over the past decade, numerous analytical solutions of 
crack problems have appeared in the open literature [1]*. A 
great number of these solutions are concerned with idealized 
crack geometries in plane or axisymmetric elasticity. How- 
ever, only a few problems involving the interaction of cracks 
with neighboring boundaries have been solved satisfactorily. 

In situations where the crack intersects a free edge or sur- 
face, the method of solution becomes much more difficult and 
requires special attention. The advent of computers has no 
doubt faciliated the numerical computation of stress distri- 
butions around cracks. Without them, many of the tedious 
calculations would not be attempted. The alternating method 
is one which intimately combines analytical results with the 
numerical calculations. 

One of the requisites for solving any crack problem is to 
handle the stress singularities at the crack tips properly. 

This involves first of all a knowledge of the correct behavior 

of the stress singularity, which is a task normally accomplished 
* 

Numbers in square brackets designate one of the References 
at the end of the Chapter. 



by analytical means. Next, it is essential to preserve this 
singular behavior of the solution in the problem either by 
isolating it away from numerical computations or by treating 
it numerically with the utmost care. Generally speaking, the 
error committed near a singular point such as a crack tip or 
border will not be confined locally but will cause errors 
elsewhere as well. The same applies to corners or points in 
the elastic solid where high stress gradients are present. 

This point will be demonstrated in the present work on the 
surface crack problem in three dimensions. 

What follows is a treatment of the alternating method as 
applied to solve edge crack problems in two-dimensions and 
surface crack problems in three-dimensions. Although the 
surface crack solution is not complete, it serves as a good 
example for illustrating the complexities and understanding 
required to treat problems of this type. The mechanics of 
the alternating method is in fact rather rudimentary. It is 
described in the work of Kantorovich and Krylov [2] who ob- 
tained the solutions to potential problems by using successive, 
iterative superposition of sequences of solutions. Their 
illustration involves two sequences of solutions, each sequence 
applying to a particular geometry. By the alternating super- 
position of the sequences, the solution for the region common 
to both geometries may be found. The method is called the 
Schwarz-Newmann Alternating Technique in their book. In this 
work it will be referred to as the alternating method, and it 
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will be applied to solve elastic crack problems. 

As an example of the alternating method, consider a simpl 
problem with no singularities. Suppose the stresses in the 
quarter plane x^O, y>_0 are to be found for the boundary con- 
ditions 

T xy (0,y) = 0 T xy (x,0) = 0 

o x (0,y) ■ 0 o y (x,0) = a ( x ) 

One sequence of solutions, Sequence A, leads from the stress 
Oy A (x, 0) = q (x ) on the half plane y>0 to the stress a x /\(°,y) 
p(y) in particular and to all other stresses in the quarter 
plane considered. The second, Sequence B, for x>^) leads from 
a xB (0,y) = p(y) to a solution yielding a yB (x,0) = q(x). The 
sequences, A and B, may be formed to yield the solution for 
the quarter plane common to both half planes. 

From A, let q^(x) = a(x), x> 0 and q^(x) = a(-x), x<0. 
For the first solution in sequence A, 

a^(x,0) = q^(x) yields a^(0,y) = p^(y), say. 

For the first solution in sequence B, 

o^(°,y) = -P^(y) yields a^(x,0) = -q^^(x), say. 

Next in sequence A, 
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a yA^ x ’°) = c l^ 1 ^ x ) yields a X A^ 0,y ) = P^(y)» sa y 
and then in sequence B, 

a xB^ (°*y) = “P^(y) yields a^jp(x,0) = -q^(x), say. 

The sequences continue alternating this way. Both sequences 
apply to the quarter plane, and by superposition, the boundary 
values on the quarter plane are 

a y ( x , 0 ) = C7yj^(x,0) + a yB^ (x,0) + CTyjp(x.O) + a yg^ (x ,0) + ••• 

a x (0,y) = a xA^ ( 0 , y ) + 0 ^( 0 , y) + o^(0,y) + a*^(0,y) + ... 
or 

0 y(x,O) = q^(x) - q^(x) + q^(x) - q^(x) + ••• 
a x (o,y) = p^^(y) - p^^(y) + p ^ " P^ ^ (y) + 

Therefore, after superposing the first n terms of each sequence, 
the boundary values are 

a y (x, 0) = a (x ) - q (n) (x) 

a x (0,y) = 0 

Additional solutions from the sequences are superposed until 
the residual stress, q^ n ^(x), on the x-axis is negligible 
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compared to the applied stress a(x). 


A different form of the alternating method was applied to 
the problem of an edge crack in a semi -i nf i ni te region by 
Irwin [3], who essentially used only the first few solutions 
of each sequence to estimate the stress intensity factor. 
Lachenbruch [4] included additional terms from each sequence 
to obtain a better estimate. Further study of the alternating 
method applied to this problem will be carried out in detail 
in Section 2. One of the two sequences of solutions involved 
there is for an infinite region containing a finite straight- 
crack subjected to arbitrary normal stresses on its faces. 

The solution is written in integral form and may be evaluated 
for any given loads. In this case, the analytical result 
gives the singularity separately, and no additional special 
numerical treatment is required. 

Recently, application of the alternating method to finite 
width strips cracked on one edge has been formulated by Broz 
[5]. The method of Section 2 could be extended to this case 
by using three sequences of solutions. Reference to such 

V 

extensions will' be made later on. However, Broz considers 
more complicated analytical problems involving only two se- 
quences. One sequence considers a row of collinear equally 
spaced cracks of equal length, each loaded by the same arbi- 
trary normal stresses on its faces. He gives the formulas 
for the stresses at the locations of the edges of the strip 
in terms of integrals which may be numerically evaluated. 
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The other sequence gives a similar expression for the stress 
produced at the location of the crack in an uncracked strip 
loaded arbitrarily on its edges. He has included very few 
numerical results in [5], but a computer program for obtaining 
additional results could be easily written. 

In a similar fashion, three-dimens ional problems may also 
be solved by the alternating method. Smith et al . [6,7,8] 
have published a series of papers containing numerical results 
for the effect of interaction of cracks with free surfaces. 

One of the geometries considered [6] was a surface crack per- 
pendicular to the surface of a half space. The crack is in 
the shape of half of a flat circular disk with the diameter 
on the surface. The other problems involved cases where the 
crack is a larger of smaller portion of the circular disk 
whose diameter parallel to the surface lies either above or 
below the surface. In [7] the crack intersects the free sur- 
face, and in [8] it is completely embedded in the half space. 

A generalization of [8] was later made by Shah and Kobayashi 
[9] in which the embedded crack is elliptical in shape with 
the minor axis perpendicular to the surface. 

As in the plane problems, one sequence of solutions for 
the three-dimensional problems above is chosen to be an in- 
finite body containing a crack. For the penny -shaped crack 
a fairly convenient solution is known for pressure on the 
crack of the form r n cosme where r and e are the usual polar 
coordinates. A less convenient solution for the elliptical 

crack is known for polynomial variations of the stress on 
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the crack. But these solutions account only partially for 
the effect of the stress singularity. The other sequence 
used for the three-dimensional problems is made up of solu- 
tions of arbitrarily loaded half spaces. It may be seen in 
Section 3 that the stress applied to the half space must be 
singular. The previous solutions do not satisfactorily 
treat this effect, which is shown in Section 3.5 to drasti- 
cally alter the numerical results. The singularity on the 
surface of the half space may be taken care of analytically, 
and then one is faced with an additional singularity in the 
other sequence. That is, the stress applied to the crack is 
singular at two points of the crack front. No analytical 
way of handling this singularity could be found at present, 
but special precaution was taken numerically. 


2- EDGE CRACK PROBLEM 

The determination of stress intensity factors for bodies 
containing cracks extending inward from a free edge has at- 
tained importance because such cracks are frequently found in 
structural members. The analytical model is usually based 
on the theory of elasticity in plane strain or generalized 
plane stress. The case of a crack perpendicular to the edge 
of a semi - i nf i ni te material (Figure 2.1) will be considered. 

Figure 2.1 


- 7 - 



Irwin [3,10] estimated the stress intensity factor, k-j , 
for the case of uniform pressure on the edge crack (or uni- 
form tension at infinity perpendicular to the crack) and ob- 
tained a value of k-| ten percent higher than that for an 
infinite plate with a crack of twice the length. He used 
two problems whose solutions are given in Sections 2.1 and 
2.2. The infinite plate with an arbitrarily pressurized 
crack is considered first. Then the half plate subjected to 
arbitrary normal stresses is superposed on the infinite 
cracked plate. Irwin used the superposition to obtain two 
simultaneous integral equations (see Section 2.3) in two 
unknown functions and solved them by an iterative method which 
is equivalent to the alternating method presented here. 

Later, Bueckner [11] reconsidered the problem and ob- 
tained a more accurate correction (13%)*. By superposing the 
two problems described above, he obtained a singular inte- 
gral equation for the deformed shape of the crack. The equa- 
tion permitted straightforward integration for obtaining the 
pressure required to produce certain prescribed deformed 
shapes. A number of different deformed shapes and associated 
pressure distributions were obtained. A linear combination of 
these pressures was fitted to the actual uniform pressure by 
collocation (equating the actual pressure to the linear 

Winne and Wundt [12] overlooked the correction of Bueckner 
when they used some of his other results on cracks in finite 
width strips (also contained in [11]). The review of Paris 
and Sih [1] points out the error in [12]. 
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combination at a number of points equal to the number of un- 
determined coefficients). Thus the deformed shape of the 
crack was determined as a known linear combination of the 
prescribed shapes. And the stress intensity factor followed 
immediately (as from equation (2.16) below). 

Koiter [13], in his analysis of the same problem, gave a 
refined value of the stress intensity factor involving a cor- 
rection of 12.15%. The computation of a definite integral was 
done numerically with such precision that the factor above is 
in error by no more than one unit in the last digit. The 
definite integral results from the application [14] of Mellin 
transforms [15] to the quarter plane problem to which the one 
under consideration reduces by symmetry. The transform leads 
to an equation of the WIENER-HOPF type which is then solved. 

Table I: k-j values for uniform tension in Figure 2.1 


Source 

kj/a/a 

Irwin [3,10] 

1 .1 

Bueckner [11] 

1.13 

Koiter [13,14] 

1 .1215 

Wigglesworth [16] 

1.122 

Lachenbruch [4] 

1.1 

Stallybrass [18] 

1 .1215 

Sneddon [19] 

1 .1215 

Appendi x 

1.1215* 


•ff 

obtained by using equation (2.50). 
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A more general study of the problem using the alternating 
method, was made by Lachenbruch [4], who considered variable 
as well as uniform loads on the crack. 

There have been other solutions [16,17,18] of the problem 
using the same technique as in [13,14]. The latest by 
Stallybrass [18] obtains the stress intensity factors for 
stress on the crack varying as a power of distance from the 
edge. He gives numerical results for integer powers, 0, 1, 

2, ••• 10, which could be combined if the stress on the crack 
is a polynomial function of that distance. Sneddon [19] also 
used the quarter plane formulation, but his technique of 
solving the problem (similar to the method of Section 2.1) 
leads to a Fredholm integral equation of the second kind. Two 
approximate methods of solving the integral equation are com- 
pared for the case of uniform stress on the crack, and both 
give very good accuracy. 

Sections 2.1 and 2.2 which follow contain the solutions 
of the two problems used in the alternating method. The 
alternating method is presented in detail in Section 2.3 for 
the case of concentrated forces acting on the crack. The 
numerical results of Section 2.4 are used to develop an inte- 
gral expression for the stress intensity factor for arbitrary 
loads on the crack. The integral is evaluated for several 
examples which are contained in the Appendix. One of the 
examples is the problem of uniform pressure on the crack which 
has been considered by many others. A comparison of this 
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result with those obtained by the various other investigators 
can be found in Table I. 


2.1 Infinite Plate With a Central Crack 


The method of integral transforms [15] may be used 
to obtain a solution of the equations of elasticity general 
enough to satisfy the boundary conditions of this problem. 
Formulas from Chapter 9 of [15] for the special case of sym- 
metry about the x- and y-axes may be written in integral form. 
For the present work, it is sufficient to note that the ex- 
pressions satisfy all of the equations of elasticity by direct 
substitution. For the case of plane strain the displacements 

u and v in the x- and y-di recti ons , respectively, are given 
x y 

by 


E 


T+v 


E 

1+v 


00 

» 

[(1 -2v)-ay]B(a)e a ^sinaxda 

0 . 
oo 

| [2 (1 -v)+ay]B(a)e a ^cosaxda 
0 


( 2 . 1 ) 


where E and v are the modulus of elasticity and Poisson's ratio. 
B(a) is a function which will be determined by the boundary 
conditions. The stresses associated with the displacement 
field of equations (2.1) are 


a x 


= -J (1 -ay)aB(a)e a ^cosaxda 

S 

= -j (l+ay)aB(a)e a ^cosaxda 
0 


( 2 . 2 ) 
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oo 


T = -y a 2 B (a)e~ ay sinaxda 
xy J 

0 

The solution given by equations (2.1) and (2.2) may 
also be used for generalized plane stress if E is replaced by 
(1 +2v) E/ ( 1 +v) 2 and v by v/(l+v). Equations (2.1) and (2.2) 
should be applied only for yj>0. For negative y, the conditions 
of symmetry about the x-axis should be used. The crack pro- 
blem of this Section is restricted to the case of loading by 
normal stresses on the crack which are symmetric about the 
y-axis. With this restriction, the problem has the symmetry 
required for application of equations (2.1) and (2.2). It 
may be noticed that the last of equations (2.2) satisfies the 

boundary condition t (x, 0) = 0 on the crack and the x-axis 

xy 

of symmetry outside of the crack. 

Figure 2.2 

The boundary conditions which determine the value of 
B(a) are 

o y (x,0) = -q (x ) , | x | <a v y (x,0) = 0, |x|>a (2.3) 

where the first of equations (2.3) comes from the imposed 
pressure, q(x), on the crack. The second of equations (2.3) 
results from symmetry about the x-axis. 
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The boundary conditions (2.3) require these special 
results from .eqs . (2. 1 ) and (2.2): 


fTv V x ’°) = 


B(a)cosaxda 


° y (x,0) = 


(2.4) 


aB (a)cosaxda 


Let v(x) = v w (x,0) = [ B(a)cosaxda 

2 ( i - v 2 ) y l 

The Fourier inversion [15] of equation (2.5) gives 

00 

B(a) = f j v(x)cosaxdx 
0 

But by the second of equations (2.3), v(x) = 0 for x>a 
Therefore, B(a) becomes 
a 


(2.5) 


B(a) = L 


v (x) cosaxdx 


( 2 . 6 ) 


Now, since v(x) is a displacement, the representation 


v(x) = 


\jj_[ t) tdt 

\. 


0<x<a 


(2.7) 


may be used. Integration of equation (2.7) by parts shows 
the correct crack tip opening shape. If equation (2.7) is 
substituted into equation (2.6), 


- 13 - 


B(o) 


a 

cosaxdx 

0 


f ^(t)tdt 


x 


/t 2 - x 2 



^(t )tdt 


f cosax_ dx 
J /t 2 -x 2 


a 

| ^(t)J Q (at)tdt 
0 


( 2 . 8 )* 


The first of equations (2.3) and the second of (2.4) 
give the equation 

oo 

| aB(a)cosaxda = q(x) , 0<x<a (2.9) 

0 


Integrating equation(2. 9) , 

00 x 

B(a)sinaxda = J q(£)d£ , 0<x<a (2.10) 

0 0 


If equation (2.8) is substituted into equation(2 . 1 0) an 

integral equation for i|j(t) results. 

00 a x 

| sinaxda J ^(t)J Q (at)tdt = j q(5)d£ , 0<x<a 

0 0 0 


a 00 

f r 


X 


or 


j i|j(t)tdt 
0 


J J Q (at)si naxda 
0 


J qU)d£ , 0<x<a 

0 


( 2 . 11 ) 


* 

The formula for J o ( x ) > the Bessel function of order zero', 
on page 27 of [20] is used in the last step. 
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But ' 


Jg(at)sinaxda 


0 

1 


^ /x 2 - t 2 


t>x 

t< X 


Therefore, equation (2.11) reduces to 


f 3^^ = f q(C)« , 


0<x<a 


( 2 . 12 ) 


Equation (2.12) may, be solved as a special case of Abel's 


integral equation 
x 


** 


( Mlldt . h(x) _ 

J Q Zx2-t2 


which has the solution 


0<x<a 


<f>(t) = | (t | hl |*^ x + h( 0 ) j , 

\ q v 1 2 ™ X 2 ✓ 


0<t<a 


Therefore, the solution of equation (2.12) is 
t 


♦<*> * S J 

0 


q ( x ) dx 

/t2-X2 


0<t<a 


(2.13) 


(2.14) 


(2.15) 


It is known that the opening displacement of the 

crack tip is given by (see equations (3.15) of [21]) 

_ 

First formula on page 37 of [20]. 

This form may be obtained from that on page 141 of [20] by 
a change of variables. 
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v y (x,0) = 2 <';' |i l k,/27 + 0 ( r ) 


(2.16) 


for plane strain. In this result, r = a-x, and is the 
Mode I stress intensity factor. From equation (2.5) 

v y (x,0) = v ( x ) 

and from equation (2.7) for 

v ( x ) = f MilMt 

a _ r /t 2 - (a-r ) 2 

In equation (2.18), when r is small, r and p may be neglected 
in comparison with a, and then 

v ( x ) = + 0 ( r ) 

/2a /r-p 

= /2r /a ^ ( a ) + 0(r) (2.19) 

Thus, by comparing equations (2.16), (2.17), (2.19), the stress 
intensity factor 


(2.17) 


x = a-r , 


a-P ^(a-p)dp 
/r-p /2a -r-p 


(2.18) 


k.j = /a i|>(a) 


is obtained. Hence, utilizing equation (2.15), 


k, = — /a 

1 7T 


q (x)dx 
y/a 2 -X 2 


( 2 . 20 ) 
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It is also necessary to evaluate the normal stress 
on the y-axis which, from equations (2.2), is 


° x ( 0 ,y ) = -[ (1 -ay)aB(a)e~ ay da 


( 2 . 21 ) 


From equations (2.8) and (2.15) 


B(o) - ? 


j (ctjtdt f 

l U ' A 2 -X 2 


0 

a 


= § j q(x)dx 
0 


J 0 (at) 

/ 1 2 - X 2 


tdt 


Therefore, interchanging the order of integration in equation 

( 2 . 21 ), 

a 

( 2 . 22 ) 


° x (0,y) = - f(x,y)q(x)dx 


where 


0 


f(x,y) = g(y»t) 


tdt 


/t 2 - x 2 


(2.23) 


where g(y»t) = 


( 1 -ay )aJg (at ) e~ ay da 


The evaluation of g(y,t) leads directly to 




(2.24) 


The integration of equation (2.23) may be accomplished with 
g(y,t) given by equation (2.24). The result is 
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(2.25) 


f (x,y) 


- y_ . , 

/ a 2 -x 2 

2y 2 2a 2 +y 2 “ 

2 2 1 

x +y 

V a 2 +y 2 

_x 2 +y 2 a 2 +y 2 _ 


Finally, substituting equation (2.25) into (2.22), the in- 
tegral 


p(y) = 


2 

71 /p+72 


a 

0 


q ( x ) 


/a 2 -x 2 F 2y 2 
x 2 +y 2 _x 2 +y 2 


dx 

a 2 +y 2 _ 


(2.26) 


for the desired stress. 


a x (0,y ) = p(y) (2.27) 

is obtained. 

The above results apply also to the case of shear 
loading on the crack. In fact if the stress on the crack is 


T xy ( x »°) = -q(x) , I x | < a 


where q(x) = q(-x), then the Mode II stress intensity factor 


i s 



q (x)dx 
/a 2 -x 2 


(2.28) 


and 


T X y(o,y) = P(y) 


where p(y) is given by equation (2.26). 



2.2 Edge - Loaded Semi - Infinite Plate 

The solution given by equations (2.1) and (2.2) 
applies to the case in which a half plane, y_>0 (Figure 2.3) 
is subjected to normal stresses 

o y (x,0) = P(x) (2.29) 

symmetric about the y-axis. It is found from the second of 
equations (2.2) and the Fourier inversion theorem that 

co 

aB(a) = - f p(x)cosaxdx (2.30) 

0 

Figure 2.3 


The first of equations (2.2) and equation (2.30) 
combine to give 


o x (0,y) 


2 

TT 


( 1 -ay )ae" a ^da 


-0 


p (x )cosaxdx 


2 

ir 


p(x)dx ( 1 -ay )ae~ a ^cosaxda 
0 


00 

4 f v 2 


= ? y J 


x 2 p (x)dx 


J 0 ( x 2 +y 2 ) 2 


(2.31) 


A change of coordinates gives us that for a half 
plane, x>0, subjected to tension on the edge. 


a x (0,y) = p (y ) 


(2.32) 
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the desired stress is 



2 . 3 I terati ve Formulation of the Problem 

The iterative approach may be though about in two 
ways which are essentially equivalent. One way used by 
Irwin [3], involves the superposition of the preceding two 
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solutions. The boundary conditions of the edge-cracked plate 
problem result in a pair of integral equations: 


q(x) - i x f r 2 P(y> d ? = o( x ) , 0<x<a 

" J 0 <x 2 +y 2 ) 2 


ana 


p(y) 


2 y 

\ 

q (x)/a 2 -x 2 

r 2y 2 2a 2 +y 2 

77 /a 2 +y 2 

) x 2 +y 2 

j_ x 2 +y 2 a 2 +y 2 _ 


dx = 0, 


0<y« 


where a(x) is the stress applied to the edge crack. These 
equations may be solved by iteration by assuming some form 
for q(x), say q(x) = a(x), and using this assumed form in the 

'■'S. 

second equation to solve for p(y). This value of p(y) is 
then used in the first equation to obtain a corrected value 
of q(x). Then a corrected value of p(y) is found, and the 
procedure is continued until there is no significant change 
in p(y) or q(x). The last value of q(x) gives the stress 
intensity factor according to equation (2.20). 

A more satisfactory way, from a physical and computational 
point of view, of interpreting the iterative approach involves 
the alternating superposition of the previously developed 
solutions. Consider the details of the case, 


a(x) = P6(x-b), 0<b<a 

where 6(x) is the Dirac Delta generalized function. The 
above expression represents a line force, P (force per unit 
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length in the z-di recti on ) , acting a distance b from the free 
edge. (Figure 2.4). 


Figure 2.4 


Let us refer to the infinite plate containing the 
crack as Problem A, and to the edge-loaded semi -i nf i ni te plate 
as Problem B. Subscripts of A or B will denote the problem 
with which the subscripted stress is associated. Then the 
iterations consist of superimposing Problem A and Problem B. 
The procedure begins with Problem A. 

In Problem A, let the stress on the crack be 

o<j}>(x,0) - -q (0) (x) = -o(x) - -pa(x-b), x>0 (2.37) 

and require it to be even in x. That is, the crack is opened 
by four symmetrically located concentrated forces. The pres- 
ent simple form of cr(x) allows 

k (0) = 2 P/F 
1 71 /a 4 -b T 

and cj^(0,y) = p^(y) 


where 


P<0)(y) * /EE [_£ iL 

77 /a 2 +y 2 ' b 2 +y 2 b 2 +y 2 


2a 2 +y 2 j 
a 2 +y 2 


(2.39) 


to be evaluated in closed form from equations (2.20) and (2.26). 
In the case of more complicated loading k j ^ and p^^(y) could 
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be evaluated by numerical integration as in the succeeding 
iterations. The second part of the zeroth integration involves 
Problem B. 


In Problem B, let the normal stress on the edge be 
o'jj'to.y) ■ -P (0) (y) 

This load gives the stress according to equation (2.33) as 
ff yB > ( x ,°) M P 0 ) (x) 

where q O)( x) = -A x ( r 2 P ( 0 ) (y)dy (2.40) 

J 0 (x 2 +y 2 ) 2 

The superposition of these two problems gives a stress state 
which satisfies all of the boundary conditions of the edge- 
cracked plate except that the stress on the crack is 

a y (x,0) = -cr ( x ) + q ^ 1 ^ ( x ) 

If q^^(x) is not negligible compared to a(x), further itera- 
tions are required. 

The first part of the first iteration requires the 
solution of Problem A for 

a^(x,o) = -q^ 1 ^(x) 

This leads to an additional contribution to the stress intensity 
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factor of 


k, (1) = I S* 


(x)dx 


and a stress on x = 0 of 


a xA^(°* y ) = p ( 1 } (y > 


where 


P ( 1 ) (y) =| 


x 2 +y 2 


: 2 +y 2 


a 2 +y 2 


An application of Problem B for a load of 


0 xB^°’ y ) = _p ^ ^ ( y ) 


gives g^(x, 0) = q^(y) 


q (2) (y) . J I d; 

71 jj (x 2 +y 2 ) 2 


The superposition of these two parts of the first 
iteration gives the solution of the edge cracked plate with 
load on the crack. 


(x ,0) = -q^^(x) + q^(x) 
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After superposing the zeroth and first interations, the stress 
on the crack is 

a y (x,0)=- a (x)+q( 2 )(x) 

and the stress intensity factor is 



The second iteration has the stresses 
ejp(x,0) = -q (2) (x) a|2 ) (o,y ) = p (2) (y) 

a x B ^ ( 0 ’ ^ ) = -P (2) (y) a yB^ x ’°) = 

f 2 ) 

and the stress intensity factor associated with q v ; (x). The 
functions p ^ 2 ^ ( y ) and q^)( x ) are given by formulas similar to 
equations(2.40) and (2.41). After superposing this interation 
on the others, we have the edge-cracked plate subjected to 
stress 

a y (x,0) = -a(x) + q ^ 3 ^ (x) 

on the crack. For this load the stress intensity factor is 



The pattern of the iterations should be clear now. 
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For the computations , introduce the dimensionless 

variabl es 


€ = x/a u = y/(a+y) c = b/a 

and the functions 

s ( n ) ( 5 ) = — 2 “ b 2 q^ n ^(a?) 

t (")(u) = /JTTBT p(n)(^_) 

Then the Iteration is from t'^'(u) to s'^'(C) to t^^(u), etc. 
where 


t (°)( u ) = 2 1 ~ c 2 u 3 (1-u) 2 |~ 2u 2 

71 c 2 (l-u) 2 +u 2 /(l-u) 2 +u 2 _c 2 (l-u) 2 +u 2 


1 


s (n) (5) = -I 5 j 

0 


t (n ~ 1) (u)du 

[C 2 (l-u) 2 +u 2 ] 2 


9 


n>l 


2 ( 1 -u ) 2 +u 2 ~1 
( 1 -u ) 2 +u 2 J 

( 2 . 42 ) 

( 2 . 43 ) 


t (n) (u) = I 


(1-u) 


/(l-u) 2 +u !i J 0 


s (n) (g)/T^ f 2u 2 

£ 2 (l-u) 2 +u 2 _£ 2 (l-u) 2 +u 2 


2 ( 1 -u ) 2 +u 2 
(l-u) 2 +u 2 _ 


dC , 


n>l 


( 2 . 44 ) 


In terms of 

s(?) = s (1 )(0 + s (2) U) + . . . 
the stress intensity factor is given by 
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(2.45) 



where 

1 

f(c) = 


The dependence of f(c) on c comes from t^(u). 



The stress intensity factor for the edge-cracked 
plate subjected to shear forces (Figure 2.4) is given in terms 
of the same function f(c) above as 


/a 2 - b* 


[1 + f(£)] 


(2.46) 


2 . 4 Numeri cal Resul ts and Green 1 s Functi on 

The numerical iteration of equations (2.42-44) is 
straightforward except when the concentrated forces are very 
near the edge of the plate. For b/a _> .2 the integrals were 
evaluated using Simpson's rule with 250 subdivisions, and five 
iterations were used. Tests using more subdivisions and it- 
erations showed that the accuracy of the results is about one 
percent. As b/a was decreased, more subdivisions were required 
and it appeared that eight iterations were required. The 
results are shown in Figure 2.5. The correction factor plotted 
there is the fractional increase due to the edge of the 
stress intensity factor for an infinite plate with four sym- 
metrically located normal or shear forces. 
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Figure 2.5 


No values of the function f(b/a) were obtained for 
b/a < .05, but the curve plotted in Figure 2.5 seems a rea- 
sonable extension of the computed points. The computed points 
are shown as boxes, and the curve drawn through the points is 

f(c) = (1 -c 2 ) [0. 2945 - 0.391 2 c 2 + 0.7685 c 4 

- 0.9942 c 6 + 0.5094 c 8 ] (2.47) 

This function will be used in a Green's function analysis of 
the stress intensity factor for an edge crack subjected to 
some arbitrary distribution of pressure. 

One reason for the difficulty in obtaining points 
on the curve for small values of b/a is the singularity at the 
point of application of the concentrated load. As long as b>0, 
the stress 0 ^( 0 , y) (equation 2.39) has no singularity, and 
its removal by the half-plane solution of Section 2.2 presents 
no difficulty. But when b=0, the stress on the y-axis has a 
singularity at the origin which requires special treatment. 

A straight application of the alternating method to this case 
would, if the numerical analysis were exactly accurate, give 
a stress on the crack, (x ,0) , with a singularity at the 
origin. Each step of the procedure would leave a stress sin- 
gularity at the origin. This points up a fundamental diffi- 
culty associated with the alternating method in more general 
cases . 
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To show the difficulty, recall that equations (2.1) 
and (2.2) give the solution for the half plane, y>0, loaded 
by normal stresses on the edge, y=0. It can be seen from 
equations (2.2) that 

tf x (x,0) = a y (x,0) 

That is, both normal stresses are the same at each point of 
the edge. Applying this result to Section 2.3 for arbitrary 
stress 

o<°>(x,0) = -q (0) (x) = -ff (x) 

on the crack, it is found that the stress, o£^(0,y) satisfies 
°(0>(°,°) = <$>(0.0) or p (0) (0) = -q <0) (0). 

Similarly, in the half plane problem, 

°1b )( 0 ' 0) ' <’xB ) (°>°> or q (1) C0) = -p (0) ( 0 ) 

And so it would continue giving 

q<°) = -p«» . ,0) . _pCl) . q (2) . .. p <3) - • • • 

where each function is evaluated at the orgin. Therefore, 
q (n) (0) - q (0) (0) = o ( 0 ) 

After the n^^ 1 iteration, the superposition of all steps gives 
the exact solution of the edge crack problem for stress on 
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the crack given by 


tfy (x ,0) = -a(x) + q^(x) 

And since 

a y (0,0) = - a ( 0 ) + a(0) = 0 

the scheme does not converge to the desired solution. That 
is, q^ n ^(x) is not negligible compared to o(x) at x=0 at least. 
But no difficulty should be expected if the stress on the 
edge crack at the edge, a(0), is zero. And in the case of a 
concentrated force as in equation (2.37) a(0) = 0 as long as 
b>0. 


In the case of shear loading, the same difficulty 
exists. The equality of the normal stresses on the edge and 
on the plane at right angles to the edge has its counterpart 
in the obvious statement about the shear stresses on the same 
planes. The same kind of details could be given for this 
case, but fundamentally the difficulty is that the shear stress 
on the crack at the intersection of the crack and edge must 
be equal to that on the edge unless the stress tensor is 
allowed to be non-symmetri c . 

Consider now the case of arbitrary stress on the 
edge crack, 

cT y (x,0) = -c(x) , 0<x<a (2.48) 
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For a concentrated force, dP = <?(x)dx, located a distance x 
from the edge, equation (2.45) gives the stress intensity 
factor 


d k i 



[1+f (x/a)] 


(2.49) 


The contribution of each part of the pressure distribution 
on the crack may be written as in equation (2.49). The total 
stress intensity factor is obtained by integrating, 



(2.50) 


For shear loading, t (x,0) = -t(x), k« is given by the same 

xy l 

expression with a(x) replaced by t(x). Note that the first 
term is the stress intensity factor for an infinite plate 
containing a crack of length 2a subjected to the same pressure 
(even in x). The second term is the correction due to the 
presence of the free edge. The results of integrating equation 

(2.50) for a number of pressure variations are contained in 
the Appendix. 

The value of the stress intensity factor listed in 
Table I was obtained by application of equation (2.50). It 
may be seen to agree with the value calculated by Koiter [13]. 
In addition, the other edge crack results in the Appendix agree 
with those of Lachenbruch [4], though he published only two 
digit stress intensity factors. The result for linearly vary- 
ing pressure agrees with that of Stallybrass [18]. 
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3. SURFACE CRACK PROBLEM 


One of the crack geometries that has received continual 
interest in fracture mechanics is that of a semi -el 1 i pti cal 
crack whose major axis lies on a stress free surface. This 
configuration is analogous to the two dimensional version of 
the edge crack problem discussed earlier except that no re- 
liable method has yet been found to solve the three-dimen- 
sional problem. The major cause of difficulty lies in the 
lack of information at the points where the crack border 
intersects with the free surface. Sih [22] has discussed this 
point in detail in connection with the finite thickness crack 
problem. This difficulty is clearly evidenced by the varia- 
tion among results [6,9,10] published on the semi-circular 
and semi -el 1 i pti cal surface crack problems. In [6,9], the 
maximum value of the stress- i ntensi ty factor was found to be 
on the free surface, whereas in [10] the maximum occurred at 
the utmost interior point of the crack front. 

In order to demonstrate the sensitivity of the solution 
to the influence of the free surface the semi-circular crack 
problem is again treated in this section by the alternating 
method. With special care given to the stress state around 
the critical points where the crack border intersects with 
the free surface, a more accurate numerical solution was 
achieved. The results are strikingly different from those 
published in the literature [6,7,8,9,10,23] and are given 
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for each iteration to illustrate the rate of convergence. 

The construction of the two sequences used in the alternating 
method is described in Section 3.3. The general solution used 
in one sequence is obtained in Section 3.1. The half space 
solution required by the other sequence is presented in Sec- 
tion 3.2. In Section 3.4, the singularities which affect 
the numerical solution are discussed, and methods of treating 
them are given. And, concluding this Section, the numerical 
results are shown. 

3.1 Penny - Shaped Crack i n an Infinite Body 

The basic solution used for the infinite body con- 
taining a penny-shaped crack (Figure 3.1) is obtained from 
Muki [24]. His solution is valid for completely arbitrary 
loading of a semi - i nf i ni te solid. When the loading is re- 
stricted to normal stresses on the surface, his solution may 
be applied to the case of a penny-shaped crack subjected to 
arbitrary variation of pressure on its faces. The crack 
problem may be reduced by symmetry to a half space problem 
in which known normal stress acts inside a circular region 
and unknown normal stress (producing known displacement) acts 
outside this region. 


Figure 3.1 
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A characteristic length is taken to be a, and all 
length coordinates are expressed as the product of a and a 
dimensionless coordinate. The cylindrical coordinates (ar, 

6, az) shown in Figure 3.1 involve the dimensionless variables 
r and z. In the following equations, a quantity, p Q , which 
has dimensions of stress is shown explicitly so that all 
quantities appearing after the summation signs are dimension- 
less. In this Section a will be the radius of the crack and 
Pq, the uniform pressure acting on it. In other applications 
they would designate some other parameters. 


If the loading is further restricted to be symmetric 
about the xz-plane (0 = 0), Mu k i 1 s solution reduces to 

oo 

2yu r = -p Q a [^cosm0|[l-2v-ez]D ni (C)e‘ ?z [J m+1 (£r)-|p- J m (£r)]d£ 


2yv Q = -p Q a f^sinm0^ [1 -2v-Cz]D m (?)e"^ z J m Ur)d£ 


(3.1) 


2yw = p n a l cosm0 
2 u m=0 


[2(l-v) + £z]D m U)e' Cz J m Ur)d£ 


0 


for the displacement components, u r » v 0 , w z . The shear modulus 
and Poisson's ratio are denoted by y and v, respectively, in 
equations (3.1). Note that the dimensions of the left hand 
sides are the same as pga. The functions D m (£) remain to be 
determined by the boundary conditions on z=0. The Bessel 
functions of the first kind, J m (£r), of all integer orders 
are required in equations (3.1). 
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The stress components corresponding to the dis- 
placement field (3.1) are given by 


o° | 

o r = P 0 l cosmej- (l-£zKD m U)e" 5z J (£r)d£ 
r u m=0 [ ' m 

oo 

+ I | 0-2v-5z)D m (C)e- 52 t!^l!- J m (5r) * J m+1 (5r)d 5 


° e = P 0 I cosme(-2vf?D U)e~^ z J m Ur)d5 
m=0 '0 


m 


- i ( (l-2v-Cz)D m (C)e- ?z [!!ii!!fi) J m Ur) + O mtl (?r)]dc 




°z = ‘ p 0 £ cosme (l+5z)£D (Oe _?z J (Cr)dC 
m=° 


( 3 . 2 ) 


W 

°o * 

T ez = " p 0 r £ m sinnie ?D m^^ e "^ Zj m^ r ^ d? 
m=0 i 


T rz = " p 0 r £ cosme ?D m (^)e‘ cz [Cr J +1 (?r) - mJ (£r)]d£ 
m=0 


T re = p 0 r L m s1nme 
m=u 


(l-2v-Cz)D m (?)e" Cz [J rn+] (Cr) 


0 


m-1 

£r 


J m (Cr)]dC 
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The shear stresses, t Qz and t^ z , vanish on the plane 
z=0 as required by symmetry about the xy-plane. The boundary 
conditions of the penny-shaped crack problem which must be 
satisfied by appropriately choosing D m (5) are 

°z ~ -P 0 P( r ’ e > on z ~® ^ or r< ^ 

(3.3) 

w z = 0 on z=0 for r>l 

where p(r,0) is the dimensionless pressure acting on the faces 
of the crack. Since p(r,0) is required to be symmetric about 
the xz-plane, it may be expanded in a cosine series as 


p(r,e) = K m (r)cosme 


m=0 


where 


K 0 (r) - 7 


0 


p(r ,0)d9 


K m (r) = f | p(r ,0)cosm0d0 , mj>l 


(3.4) 


(3.5) 


The equations for D m (£) are obtained by using 
equations (3.3) in conjunction with (3.1), (3.2), (3.4). 
They are 


H< f >V« r > dE - K r, (r) > r<1 

0 

CO 

j D m (C)J m (Cr)d? = 0 , r>l 

0 


(3.6) 


for each m = 0, 1, 2, 


• • • 


The procedure for solving these 
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equations is analagous to that in Section 2.1. The displace- 
ment-like function, 


w ( r ) 
m ' 


D m U)J m Ur)d£ 
J m m 

0 


(3.7) 


is introduced. Since by the second of equations (3.6), 
w m (r) = 0 for r>l , the Hankel inverse transform [15] of 
equation (3.7) is 
1 

D m U) = ? | w m ( r ) J m (S r ) dr (3.8) 

0 


The displacement-like behavior of w m (£) is assured if it is 
represented in the form 


w m ( r ) 
m 


i r- 

7 T 


9m (t > 


/t2-r2 


(3.9) 


Use of equation (3.9) in equation . (3 .8) gives another repre- 
sentation of D (£) , 




h 


J m (Cr) r m dr 


0 

1 


9 m U) 


'm 


t-^dt 

7£TT7T 


= I c f t" m+1 g (t)dt | r m+1 Jm(€r)dr 

77 l 01 l /r—f? 


i 

= - £ [ t 2 g m (t) j m (£t)dt 
n 3 m m ^ 


. . * 
(3.10) 


* ^ 
The formula at the top of page 32 of [20] may be used in 

the last step if the substitution, r=tsin0 is made. 
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where 


V 5t > ■ W 5t) 


is the spherical Bessel function of the first kind of order m. 


If the first of equations (3.6) is multiplied on 

m i I 

both sides by r m , it can be integrated with respect to r 
to become 


oo 

r 


D m <e)r- m+1 J m+1 Cer)d ? - L m (r), 


r<l 


r 

where L m (r) = j n m+1 K m (n )dn 

0 


(3.11) 

(3.12) 


The equation' which determines g m (t) is obtained by substituting 
equation (3.10) into (3.11): 


2 

TT 


S0 m+1 (5r)dC 


0 

1 

I 

0 


t 2 g m (t)j (£t)dt = r” m "H (r), r<l 
3 m v 7 J m v m ' 


0 


(3.13) 


or j t2 9 m ( t ) dt ?J m+l ^ r ) J m^ t ) d? = \ r " m ' lL m( r )’ r<1 


Knowing that 1 


oo 

r 


?J m + l ( ^ )j m Ut)d ^ = < 


0 


t m r -m-l 


, t>r 
, t<r 


^ /r 2 - t 2 

equation (3.13) may be written as 

/ + m+2 .. 

I 9 m (t) — - | L <r) , 

1 m /Fr-I* ^ m 


r<l 


(3.14) 


Third formula on page 37 of [20] with y=m, v=m+^ 
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By application of the solution, equation (2.14), of equation 
(2.13) to equation (3,14), 


t , « 


t m+2 g m (t) = t 


L (r)dr 
m 

✓t 2 - r 2 


or, since equation (3.12) gives 


L_( r ) 
m 


r m+1 K m (r) 
m ' 




= t 


' m - 1 f 
: 0 


r m+1 dr 
/t 2 - r 2 


(3.15) 


Changing the variable, of integration in equation (3.15) from 
r to t, = r/t gives the form, 


«.<*> * I 


m+1 

/T^EJ 


The stress intensity factor is most conveniently 
determined by examining the crack opening shape near r=l . 
From equations (3.1), (3.7) the displacement for 2=0 is 


w_ = 


1-v 


z , y 


P 0 a 1 w m (r)cosm0 
m= 0 


(3.16) 


For r=l-p, equation (3.9) becomes 


w m (l-p) 


f o-p) 


m 


+ -m+l j . 
*.<*> L - 


i-p 

P 


/t 2 -r 2 


* I | V 1 '"’ 


/p-n /2-p-n 
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As p+O, p and n may be neglected in comparison with 1, and 
this leads to 


w ( 1 - p ) = f g (1 ) -1 f -S3- 

m * " /? i /Fn 


+ 0(p ) 


= | 9 m O + °^ p ) 


Therefore, the equation (3.16) for the crack opening near the 
crack front becomes 


W 2 = f ^ 4 — Vo/a /2a? I g (V)cosme + Ofo ) (3.17) 

m= 0 

The relationship between displacement and stress intensity 
factor (See equation (5.37) of [21]) is the same as in equa- 
tion (2.16) for plane strain. In equation (2.16), note that 
r is not dimensionless. Hence it is identified with ap in 
equation (3.17). Thus, comparing equations (3.17) and (2.16), 

k = £ P 0 ^ I g (1 )cosm6 (3.18) 

P tt u m=0 m 

In equation (3.18), kp is the physically dimensioned stress 
intensity factor. The subscript p is used to distinguish it 
from the dimensionless stress intensity factor 

00 

k ( 0 } = I g m (l)cosme = k /(£ p Q /a) (3.19) 

m=0 H u 

which will be used in the numerical calculations. 
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The stress intensity factor may also be determined by 
calculating the singular part of the stress distribution. If 
equation (3.10) is integrated by parts, 


D (g) = - 
nr ' 77 


^m^ ^ J m+1 ^ ^ 


2 

7T 




m 


g m (t)]dt 


only the first term contributes to the stress singularities. 
The spherical Bessel functions may be written as the product 
of polynomials in ^ and trigonometric functions as e.g. 


j 0 u> - 


si ng 
g :■* 


^(O 


sing 

g 2 


cosg 

g 


The resulting integrals of the type 

oo 

p 

. g n e" Cz J m (gr)singdg 
0 

obtained from equations (3.2) may be evaluated from formulas 
based on the third one on page 33 of [20]. By differentiating 
and separating the real and imaginary parts of the formula 
referred to, results of the type 

oo 

f e' ?z J (gr)singdg = ^ sin(f-+ f^) + 0(1) 

J 0 m /Zp 2 2 

are found. After these results are used, and the nonsingular 
terms discarded, the stress state may be written in the form. 


| P 0 Ccosf - ^ sin4>sin^|-] + 0(1) 


/Tp 
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(3,20) 


2 _ k(l) 

/2p 


= 2V ? p 0 ~ [cosf] + o(i) 


= § P n ~~ [cosf + sin<f>sin^|] + 0(1) 


/Zp 


T rz = 'I p 0 ^ sin<f>cos^4-] + 0(1) 


/Fp 


In equations (3.20), p and <p are the dimensionless polar 
coordinates shown in Figure 3.2. The stresses x Qz and x Q , 
are nonsingular. It may be verified that the singular terms 
satisfy the condition of plane strain, 


°e = v(o r + °z> 
locally. 


Figure 3.2 

The present iterative aprroach requires the calcula- 
tion of the stresses on the yz-plane. In general each of 

the stresses, a , x , t , is not zero, but for the case 

x xy x l . 

in which the pressure on the crack is symmetric about the 
yz-plane, only the normal stress is not zero. It may be 
evaluated from equations (3.2) by 

a x (0>y>z) = o 0 (y,7T/2, z) 

It is helpful in the numerical analysis to separate the 
singular and nonsingular parts of this stress as 
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a 

x 




] + p 0 q(y>z) 


(3.21) 


where (Pj*^ ) and ( P 2 » 2 ) are coordinates in the yz-plane 
as shown in Figure 3.4. The nonsingular function q(y,z) is 
computed from the preceding two equations as long as the 
point ( 0,y ,z ) is not close to the crack front. When (0 ,y , z ) 
nears one of the points (0, + a , 0), the singular parts of 
equations (3.21) and (3.22) are equal, and q(y,z) may be 
written in a form in which the singularities are analytically 
cancelled. From equations (3.2) and (3.21), 


q(y.z) 


- k(f)[JL_ cos 






cos 



+ 




( 1 -2v-£z ) D 2 m 


C D 2 iT,(C)e-« z J 2 m (Cy)d 5 

( ? , e -?Z[M|azD. J 2m ( 5y) 


0 

+ J 2 m + 1 (5y)]de\ (3.22) 

Use has been made of the fact that for symmetry of loading 
about the yz-plane, all K, g, D with odd subscripts vanish. 


To illustrate the lack of singularity in q(y, z), the 
form valid for z = 0 and y<_l will be given. A similar expres- 
sion can be written for y^> 1 , and at y = l , both give the 
same value of q(y,z). The Fourier coefficients of p(r,6) can 
be written as polynomials, 
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, mM 


(3.23) 


K o (r) = p c + J 


K 2 m (r) - l_ A 


nil m+1 ' n 


where p c = Kq( 0) = p(0,e) is the dimensionless pressure at 
the center of the crack. For K 2m (r) in this form, equation 
(3.15) gives 


g 0 (t) - p c + J l c ln t n 

- J, Wn*" ■ (3 - 24) 


where C 

mn 


^m+n-l A mn 


(3.25) 


and I Q = J , 


= 1 


k-1 

k 


k-2 


(3.26) 


The part of Dq(£) resulting from the constant part of gg(t) 
may be evaluated as 


2 

IT 


K 


0 


t 2 p c J 0 Ut)dt = 


P c j,(C) 


and the integrals of this part of D Q (C) in equation (3.22) can 
be determined in closed form. The remaining part of D Q (£) 
and the other D^U) can be found in terms of spherical Bessel 
functions and the sine integral function. In general, the 
remaining integrals in equati on (3 . 22 ). are evaluated numerically, 
but when z=0 they can be explicitly obtained. 
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For z = 0 and y<H > the integrations in equation (3,22) 
result in 

q(y, 0 ) = I P c F 0 (y, 0 ) ♦ |j ) (-') m j 1 i 2 m+n c mn y n D- 2 v 


+ 2v ( 2m+n ) + (1 ~2v) ( 2m-2 ) ( 2m-3 ) F mn (y ) ] (3.27) 


where 


f ‘ t 4m-4 2m-4+n 

* ▼ y 


F mn^ 


= < 


I 


2m+n 


2m-4-n 


ln 2 + + TT + 


n f 2m-4 


+ (4m- 5) (4m- 4) ’ n “ 2 " 1 ’' 4 


A similar but more complicated result can be found for y>l . 
Using the polar coordinates shown in Figure 3.4, the other 
function in equation (3.27) may be written 


1 . ^ a + 4> 


F g (y,z) = 2v[ ■ sin 


j P . 


O — - — — cos - — — cos 

<- rrr ~ — c r*s c. 


/Tp. 


/2p" 


1 +2v . , 2 x A 1 - 2v / r • ^i + ^ 2 

2 'P +P o ..2 i 2 L 2 

1 2 


+ z COS 


+ --- 


2y 


+ (p^-Ucos 


y ^Pjpj 
< l > 1 +<}>2 

9 J 


[ z sin 


4* i +< J ) : 


(3.28) 


For z=0, equation (3.28) reduces to 


F 0 (y.o) = - 


v /2 ( y - 1 ) 

^2 (y+1 )+y+l 


1 + 2 V . 11 -> 2 V r—x — ■- 

T ~ 2 — arcsin - + ^T-l 


(3.29) 
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for y>J . For y<l , FQ(y,0) is equal to the constant obtained 
from equation (3.29) by putting y=l . 


Reiterating the important results, if p(r,6) is 
symmetric about both the xz- and the yz-planes, then it can 
be written 

co 

p(r>e) = - •— a (x,y,0) = J K ? (r)cos2me (3.30) 

p 0 z m=0 


It follows that 


92m (t) = 


1 

r ? 2m+l 
0 


K 2m^ t)dt 


(3,31 ) 


can be computed to yield the stress intensity factor 


oo 

Me) = f PO^ l g 2 m^ ) cos2m0 

P tt u m=0 zm 


and 


D 


2m 




t2g 2m( t)j 2m^ t ) dt * 


(3.32) 


(3.33) 


Once D 2m (£) is computed, any of the stresses and displacements 
may be determined. In particular a (0,y,z) of equation (3.21) 

A 

is obtained by calculating the function q(y,z) from eq.(3.22). 


3 . 2 Surface Loads on Half Space 
( 1 ) Nonsi nqul ar stress 

Consider the case of normal stresses a acting 

/N 

on the surface occupying the yz-plane as shown in Figure 3,3. 
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For the case 


x ( 0 ,y , z ) = p Q q(y,z) (3.34) 

in which there is no stress singularity, there are a number 
of possible methods of solution. For example, the solution 
at the beginning of Section 3.1 could be applied or the solu- 
tion for a concentrated force on a half space could be used 
as a Green's function. For the convenience of numerical cal- 
culation as in [6], the surface will be divided into squares 
sufficiently small in size so that the stress is nearly uni- 
form on each. Next a solution provided by Love [25] is used 
to obtain the contribution of the stress o on each square 

A 

to the stress o z on the xy plane. 

Figure 3.3 

Love's solution gives the normal stress a z on 
the xy-plane due to a unit tensile stress on a square of side 
b centered at (0,y,z) as 

<J z (x,y,0) = F( | y-y I > |z| , | x | ) 

where 2 tt F ( ^ , n , C ) = ( 1 -2v ) F i ( ? ,n , C ) -4 F^C.mc) 

- 2vF 3 U,n,c) (3.35) 

and the functions F.(j=l,2,3) stand for 

J 
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1.,’ 1 H + - 1.-' 1 


— & ♦ tan-' |i| - tan'’ ^ . tin -l ^ 


b-n 


[kl + b±i] + _kt!L 


2 ( b-n ) 2 +? z L a i b 2 J (b+n ) 2 +? 2 L d ‘ Cs 


j-b^i + b+i] ( 3 . 36 ) 


= tan 


-1 


Ca 


(b-^) (b-n) 


cc 


+ tan 


-1 


^ b 2 

( b+"4) ( b”nT 

Cd 


+ tan (b+5)(b+n) + tan (b-€) (b+nT ' i 


2tt , | c | <b , | n | <b 
0, otherwise 


with the following contractions 


a i = (b-U 2 + (b-n) 2 + ? 2 
b 2 = (b+£) 2 + (b-n) 2 + S 2 
c 2 = (b+5) 2 + (b+n) 2 + s 2 
d 2 = (b-£) 2 + (b+n) 2 + ? 2 


By taking advantage of symmetry, attention may 
be confined to the centers (0,y.,z.) of the squares which lie 

1 J 

in the first quadrant of the yz-plane. The effect of the 
other three symmetrically located quadrants is included in 
the result 

a z (x ,y ,0) = p Q t(x,y) 

where 

t(x.y) = I q(y 1 .i j >[F(|y 1 -yM 5 jl = M> 

i=l j=l J J 

+ F ( 15,+y | » I 5 j I . 1 X I ) ] (3.37) 
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The symmetry of the stress on the yz-plane insures that there 
will be no shear stresses on the xy-plane. 

Figure 3.4 

( 2 ) Singular stress 

When the stress on the yz-plane is singular, it 
is desirable to seek a closed form solution. For the case of 
interest, from equation (3.21), 

d\) "irl $ , I ^9 

O (0,y,z) = ^ P 0 k(J)[— cos 2 ~ + cos «“] (3.38) 

/2p i /2p 2 

such a solution can be found by applying Muki's half-space 
sol ution, equati ons (3.1) and (3.2) } to each term. The result 
i s 

■J.(x,y,0) =■ ^ p 0 k(i) [— ! S(* ) + — S(*.)] (3.39) 

where the dimensionless coordinates, p, <J>, R, iJj are shown in 
Figure 3.4. The function S ( ) is in the form of an infinite 
series involving hypergeometric functions. Rather than give 
the formulas, the results are shown graphically in Figure 3.5 
An excellent fit of the exactly computed points of the curve 
is obtained by the approximation, 

S(\p) = /$ [0.344 + 1.851^ - 1 1.637^ + 34 . 54 5ip 3 

- 61.286^“ + 68 . 1 24i|) 5 - 47 . 644^ 6 + 20.32H 7 

- 4 . 823^ e + 0 . 488ip 9 ] (3.40) 
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This expression is valid only for 0<ip<Tt/2 and for Poisson's 
ratio, v = 0.30. 


Figure 3.5 

It should be noted that there are singulari- 
ties in a (x,y,0) at the points where the crack front inter- 
sects the free surface. The numerical difficulties due to 
the singularity will be discussed in more detail in Section 
3.4. 


3 . 3 Iterati ve Formulation 

The solutions of the two preceding problems will 
be superposed alternately to solve the case of a uniformly 
pressurized half penny-shaped crack perpendicul ar to a free 
surface (Figure 3.6). The procedure begins with the well 
known solution [26] for uniform pressure of a penny-shaped 
crack in an infinite medium. This solution gives stresses on 
the yz-plane which are removed by superposing the half space 
solution for equal and opposite stresses on the yz-plane. 

As a result of the half space solution, there is a residual 
normal stress on the crack in addition to the uniform pres- 
sure. This is eliminated by superposing the solution of 
Section 3.1 for pressure on the crack equal to the residual 
stress. Again, stresses on the yz-plane are produced which 
are removed by application of the half space solution. The 
cycle of superpositions is continued until the resudual stress 
on the crack is negligible compared to the uniform pressure. 
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Figure 3.6 


As outlined, the first step in the procedure re- 
quires the solution of the penny-shaped crack problem for the 
case 


f (0) 

zA 


(x ,y , 0 ) = -p Q 


(3.41) 


Either [26] or Section 3.1 may be used to find the stress 
i ntens i ty factor , 


k p ° ^ = | P 0 ^a or k ^ 0 ^ ( 0 ) = 1 


(3.42) 


and the stress on the yz-plane, 


a xA^ ( 0 >y » z ) = P 0 k ^(J)[— — cos + — — cos tr-J 


/2p" 


/2p 


+ p 0 q^(y>z) 


(3.43) 


See Figure 3.4 for the coordinates, p,cf>. For the case of 
uniform pressure, q^(y,z) is equal to ^ FQ(y,z) which is 
given in equation (3.28). 


The second half of the zeroth iteration is used to 
eliminate the stress on the free surface. The stress applied 
to the surface is 

°xB> (0>y ’ 2 ) ■ ff^ho.y.z) (3.44) 


which gives the result (See Figure 3.4 for R,^.) 


y (0) 

zB 


- 4^ n . (0),iu r ^\l , S( ^ , 2 ) 


( x ,y , 0 ) = f- p 0 k'^(^)[ 


/nr 

1 


/nr 


] + p Q t^ ^ ( x ,y ) 
(3.45) 
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where t^(x,y) is determined by applying the method of 
Section 3.2(1). It depends only on q^, and not on the 
si ngul ar part of . 

After the second half of this iteration is sub- 
tracted from the first half, the surface will be free of 
stress, but the stress on the crack will be 

o z (x,y,0) = -p 0 - a zB ^ ( x »y 5 0) (3.46) 

rather than constant. An additional iteration is used to 
eliminate the stress 

For the first iteration, the stre.ss on the penny- 
shaped crack is taken to be 

) S(* ) 

- + 3 _] 

/2TT 

2 

+ Pg t ^ 1 ^ ( x »y) (3.47) 

Because of the singularity in the first term of equation 
(3.47), care must be taken in the numerical procedure. The 
description of the numerical method will be given in the 
following section. The contribution of the two parts of 
a zA^ are com P u ted separately and the desired quantities ex- 
pressed as 

k (1) (e) = k<°)(f) ks ( 8 ) + k* 1J (e) (3.48) 

q (1) (y,z) > I s k (0> (f)q s (y,z) + q (1) (y,z) (3.49) 

71 c b N 


r (D 

zA 


(x ,y ,0) = o 


( 0 ) 

zB 


(x.y.O) - ^ 


P 0 k ( 0 ) (f)[ 


$(♦, 
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1 


o'VtO.y.z) = p 0 k (1 ><f) [— ! — cos p- + -J — cos £*-] 


+ p n q (y» z ) 


/*p7 


/?pT 


(3.50) 


where k s (0) and p Q q s (y,z) are the dimensionless stress inten- 
sity factor and nonsingular part of o (0,y ,z) , respectively. 


due to stress on the crack of 


<*,(x,y ,0) = P n C- 


S(^j) S(i^ ) 


/2^7 




(3.51 ) 


Similarly, ki '(0) and p n q). '(y,z) stem from the stress 


p n t^(x,y) on the crack. 


The effects of these two terms are kept separate 
into the second half of the iteration. The stress on the 
half space is taken to be 


°xb > ( o>y ,z > ■ c'i'fo.y.*) 


This results in the residual stress on the crack, 


(x,y ,0) 


^ Pok ( n (|)[ f^ + fiM ] + Pnt <2) (Xi: 


/2R7 


(3.52) 


where 


t (2) (x,y) = P- k (0, (f)t s (x.y) + t< 2) (x,y) (3.53) 


In a way similar to the definitions of q $ and q^ , t $ and * 
t f. ^ are defined to be the stress a (x ,y ,0) due to the stress 

M ^ 
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a x ( 0 ,y , z ) equal to q $ and q^, respectively, and are found 
by the method of Section 3.2(1). 


Again, the second half of the iteration when sub- 
tracted from the first half frees the surface of stresses. 
Adding the first iteration to the zeroth cancells the re- 
sidual stress cr^g^ and gives 

a z (x,y,0) = -p Q - o^(x,y,0) 

Further iterations are added until the stress on 
the crack is reduced to the applied stress -pg to within 
some arbitrary tolerance. The additional computations which 
must be made for each iteration are straightforward. For 
the second iteration, the dimensionless intensity factor, 
k^ 2 ^(e), and stress a (0,y,z) = p Q q^ 2 ^(y,z) due to stress on 
the crack of a (x,,y,0) = Pgt^ 2 ^(x,y) are obtained. Then the 
half space solution gives a (x,y,0) = Pgt^ '(x,y) due to 
a x ( 0 ,y , z ) = Pgq^ 2 ^(y,z). The contribution of this iteration 
to the dimens ionless stress intensity factor is 

k (2) (e) = k (1) (f)k s (e) + k' 2) (e) 

and the stress on the crack is 


<J z (x,y,0) = -P Q - cr^B ^ ( x ,y ’ 0 ) 


where 


b (2) 

°zB 


S(4>) S(* 2 ) 

L_ + 

'mr /tft 


^ P 0 k (2) (J)[ !- + —=-1 + P 0 t' J '(x,y) 


+ ( 3 ) 
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t^(x,y) = ^ k ^ 1 ^(J)t s (x,y) + t^ 3 ^(x,y) 


For the n— iteration, the functions, k<.(0), 

P 5 (y,z), and t<j(x,y), are used again in conjunction with 

the results of the previous iteration, k^ n "^(e) and 
/ \ 

t' '(x,y). The method of Section 3.1 gives the new functions, 
k^ n )( 0 ) and q^ 0 ^ (y , z ) , which are, respectively, the dimen- 
sionless stress intensity factor and the nonsingular part 
of o ( 0 ,y,z) produced by application of the stress o _(x,y, 0 )= 
t^ n ^(x,y) to the crack. The half space solution of Section 
3.2(1) gives t^ n + ^ (x ,y ) , the stress c z (x,y,0) produced by 
a x (0,y,z) = (y,z) acting on the surface. The new addi- 

tion to the dimensionless stress intensity factor is 


k (n) ( 0 ) = p- k (n_ 1 ) (^)k s ( 0 ) + k< n) ( 0 ) 


IT 


and the residual stress on the crack is 


o'S’fx.y.O) - p P 0 k ( ">(f)[^- ♦ ^3 ♦ P 0 t( n+, >(x. 


y / 2R 1 /2RF 


.y) 


where the term carried into the next iteration, if necessary, 
i s 

t (n+1) (x,y) = p- k< n -')(f)t s (x.y) + t< ntl) (x,y) 


N 


The alternating method applied in this Section is 

modified by the separate treatment of the singular part of 

the stress (equation (3.43)) applied to the half space. The 

solution for the half space subjected to the singular stress 
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is given by the first part of equation (3.45). The 

resulting stress which must be removed from the crack, eq. 
(3.47), has two parts which are kept separate. The singular 
part is the same in each iteration except for the coefficient 
involved. Therefore, the numerical solution for a penny- 
shaped crack subjected to the singular stress (tensile) of 
equation (3.51) is determined to as high a degree of accuracy 
as possible. The contribution of this solution to each 
iteration is contained in the functions with a subscript, S. 
The remaining contribution, from the rest of equation (3.47) 
gives the terms denoted by the subscript, N. The parts are 
combined as in equations (3.48) to (3.50) to give the total 
solution of the penny-shaped crack problem for this iteration. 
The solution is treated by the same method in each successive 
iteration. 


The next section contains a description of the 
numerical details and the difficulties associated with the 
stress singularity of equation (3.47). 

3 . 4 Numer i cal T reatment of Singularities 
( 1 ) Half space 

For the case of uniform pressure on the crack, 
the first step requiring numerical analysis is the clearing 
of the stress, eq.(3.43) from the surface in the second half 
of the zeroth iteration. For convenience the problem is a 
half space with tensile stresses, 

x ( 0 ,y , z ) = p Q q(y,z) 
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on the surface. The function q(y,z) can be any of the non- 
singular functions of the first line of Table II. The solu- 
tion for the singular part of the stress a (0,y,z) is given 

X 

analytically by equations (3.38) to (3.40). 

The surface of the half space is divided into 
a number of regions as shown in Figure 3.7. Each region is 
then subdivided into a number of squares, larger squares 
being used where the stress becomes small and less rapidly 
varying. The values of q^(yjz) at the centers of each of 
the squares are computed from the closed form expression, 
equation (3.28). Because of the lengthy calculations needed 
to obtain the values of the other q's, they were computed 
at fewer points and the six point bivariate interpolation 
formula was used to get the values at the remaining points. 
The numbers in parentheses in Figure 3.7 give the number of 
values of y and the number of values of z, respectively, used 
to obtain the grid points at which q was calculated exactly. 
After the values of q(y,z) have been computed at the center 
of each square, the solution in Section 3.2(1) is used 
directly to obtain the values of the stress on the xy-plane 
denoted by 

z (x ,y , 0 ) = p Q t(x,y) (3.54) 

where the t's corresponding to the various q's are shown in 
Table II. 
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Table II: Results of Half Space Solution 

of Section 3.2(1). 



n>l 

surface stress 
a x (0,y,z)/p Q = 

q ( 0 ) (y,z) 

q s (y >z) 

q^ n) (y’ 2 ) 

results in 
o z (x,y,0)/p Q = 

t (1) (x,y) 

t $ (x,y) 

tj, n + 1 ^ (x,y) 


The values of t(x,y) are computed at the 36 
points (.2m, .2n), m , n , =0 , 1 , • • • 5 . The six point bivariate 
interpolation formula is used to obtain intermediate values 
of t(x,y). This completes the description of the half 
space problem's numerical solution. 

( 2 ) Penny-shaped crack--nonsi ngul ar load 

The solution of the penny-shaped crack problem 
is needed for stresses on the crack which contain no singu- 
larities. 


° 2 (x ,y ,0) = p Q t(x,y) 

The function t(x,y) may represent any of the t's in Table III. 
In a case of non-uniform pressure on the crack, it could also 
represent the initial stress on the crack, t^^(x,y). (In 
this example, t^(x,y) = -1). As indicated above, t(x,y) 
is known at 36 points and interpolation gives its value at 
any other point on the crack. 
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The solution of Section 3.1 requires first the 
computation of the Fourier coefficients of 

p(r,e) = - 1 (x ,y ) 
which are 

tt/2 

K 0 (r) = | | p(r,0)d0 
0 

ir/2 

K 2m (r) = ^ | p(r,0)cos2m0d0, mH (3.55) 

0 

K 2m+1 ^ r ) = 0 » m =° 

The first 40 K's are each computed for 30 values of r using 
a technique given by Filon [27]. His technique gives as 
good accuracy with 100 divisions of the range of integration 
as the usual Simpson's rule does with 500. 

The coefficients of a five term polynomial 
approximation to each K are determined by the method of 
least squares. This gives 

V r > 'Pc + j, A ln r " 

K 2n.< r ’ = n I 1 A m +1 .n' r ' n ’ < 3 - 56 > 

where A mn , m=l , 2 , • • • 40 ; n = l,2,***5 are known constants and 
P c = p(O,0). The polynomial form of the K's enable 9 m .(t) to be 
determined also in polynomial form as in equation (3.24), 
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(3.57) 


9 0 (t) = P c + j, C ln‘ n 

92m<t) = n-1 C| " +, - ntn ' 

where C mn is given by equations (3.25) and (3.26). 

The stress intensity factor may be obtained 
simply from equations (3.19) and (3.57). The integrals of 
equation (3.10) which give D m (£) may be evaluated in closed 
form in terms of spherical Bessel functions and the sine 
integral function. Recurrence relations between the inte- 
grals are used for the numerical calculations. The details 
are lengthy but straightforward. Finally, the integrals 
in eq.(3.22) involving D m (£) must be evaluated to obtain the 
function q(y,z) of equation (3.22) which gives the nonsingu- 
lar part of a (0 ,x,z)/Pq. These integrals are evaluated by 
Simpson's rule using intervals of length 0.25 and integrating 
over segments of length 5.0 until the integral over the last 
segment is small compared to the total integral. As pointed 
out in the first part of this section, q(y,z) is computed 
for various points of the yz-plane. The notation used for 
the results of this part of the numerical analysis is sum- 
marized i n Tabl e I I I . 
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Table III: Stress Intensity Factors and Stresses 

on yz-Plane for Various Loads on Crack, 



n>_l 

stress on crack 
CT Z (x ,y , 0)/p Q = 

eq. (3.51 )/p Q 

t (0) (x,y) = -1 

t 1 " 1 ( x ,y ) 

produces stress 
tf x (0,y ,z)/p Q 
whose nonsing, 
part is = 

q s (y ,z) 

q^(y.z) = 

1 F o (y » 2) 

(see eq . ( 3 . 28 ) 

qj n ^(y.z) 

produces stress 
intensity factor 

y e >4 Po 71 ° 

k s (e) 

II 

CD 

O 

(e) 

For n>2, t (n) (x,y) = ^ k (n_2 } ( |)t $ (x ,y ) + t^ n) (x ,y ) (See T. II) 
For n>l, k (n) (9) = ^ k (n_1) (J)k s (e) + kj n) (6) 


( 3 ) Penny-shaped crack -- singular load 

The singular stress on the crack was given in 
equation (3.39). For the numerical analysis consider the 
dimensionless form, 


a z (x,y ,0) 


S(*.) S(ip ) 
L_ + L_ 

/Z1T m T~ 

1 2 


Note that since o z > 0, the singular part of the load on the • 
crack gives a negative stress intensity factor. Further, 
since larger tensile stresses occur near’the intersection of 
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the crack front with the surface, the stress intensity factor 
should be expected to be larger in magnitude there. Thus, 
when combined with the other contributions to the stress in- 
tensity factor, this will cause a decrease in the total as 
the surface is approached. 

Because of the singularity, it would be advantageous 
, to have a closed form solution for this part of the problem. 
But since none could be found, it was decided to use the same 
numerical method given for nonsingular loads. It is evident 
that more terms and finer mesh sizes are required. Limita- 
tions on computer time forced the use of the same number of 
points in r and the same number of terms in the r expansion 
as above. More terms were included in the Fourier series, 
and smaller subdivisions were used in the integration. It 
was found that K m (l) decreased very slowly as m increased. 

For m=400, K m (l) was of the order of 10"^, K m (.967), of 10"^, 
K m (.933), of 10~ 4 , and for r<.5, K m (r) was of the order 10"^. 
For m=1600 these figures decreased by a factor of 10. Each 
time additional terms were included in the Fourier series, 
the magnitude of k s (e), the dimensionless stress intensity 
factor due to the singular pressure on the crack, increased 
near 0 = tt/2 and remained the same elsewhere. Finally k^(e) 
was calculated using 841 terms in the cosine series. 

The function q^(y,z) which gives the nonsingular part 
of the stress a (0,y,z) due to the load on the crack 

A 
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was computed using 241 terms in the Fourier series. This 
part of the numerical solution required more than half of 
the computer time used. But note that the same function 
is used again in each iteration. 

In reference [6] which also treats this pro- 
blem there is no separation of singular and nonsingular 
terms in the numerical analysis. The half space problem 
utilizes stresses at points within Regions I and III of 
Figure 3.6. Outside these regions, the stress is approxi- 
mated by zero. In [6], the stress is computed at 80 points 
within Region I, as compared with 45 in this analysis. The 
additional points should not be expected to account accurately 
for the singularities. 

To obtain the stresses in the regions mentioned 
above and to obtain the stress intensity factor for a given 
stress on the penny shaped crack in [6], the stress on the 
crack is computed for five values of r and 19 values of 6. 
Weddle's rule is applied to the integrals for K m (r) for 
m=0 , 2 , 4 , • • • 1 0 . Simpson's rule using, say, 40 values of 0 
should be as accurate. K m (r) 1 s obtained in [6] as a fifth 
order polynomial in r by forcing the polynomial to pass 
through the computed points. In this work, K m (r) is known 
at 30 values of r, and the fifth.-order polynomial is obtained 
by least squares curve fitting. 


- 63 - 



3.5 Discussion of Numeri cal Results 


The result. of the analysis described above is the 
curve in Figure 3.8. which shows the variation of the stress 
intensity factor along the crack front. The deepest point of 
the crack is 6=0. Initially, moving away from this point, 
there is a gradual increase in stress intensity factor. Then 
about 20° from the surface there is a rapid increase. Fol- 
lowing this a more rapid decrease occurs in the last 3°. 

Note that the stress , i ntensi ty factor tends to zero at the 
surface, a result discussed in several previous papers [22, 
28,29]. It appears that an increase in numerical accuracy 
would further decrease the value at the surface toward zero. 
However, the lack of an exact solution of the problem of a 
crack loaded by the stress of equation (3.51) is reflected 
by an accumulation of numerical inaccuracies as the number 
of iterations is increased. 

Figure .3.8 

It should be clear that the results in Figure 3.8 
are by no means the exact solution of the problem but are 
presented as the best approximation available. The trend of 
the curves differs substantially from those obtained in [6] 
using basically the same numerical procedure. The results of 
[6], shown as a dashed line in Figure 3.8, would be expected 
to differ slightly from the present results because they were 


obtained for a Poisson's ratio of 0.25 and the present ones, 
0.30. The apparent reason for the qualitative difference 
between the results is the differing treatments of the sin- 
gularities. As background information some intermediate 
results will be given. 

The curve of Figure 3.8 is the result after the 
fourth iteration. The stress intensity factors after each 
iteration are shown in Figure 3.9. Each new iteration raises 
the stress intensity factor, but the amount of increase 
decreases as the number of iterations increases, with one 
exception. Near the surface the curve for the second itera- 
tion lies below that for the first. 

As described in Section 3.3, the contribution of 
each iteration to the stress intensity factor consists of 
two parts. The two contributions for the first iteration are 
shown in Figure 3.10. The function k s ( e ) is used with each 
iteration as the stress intensity factor for the singular 
pressure on the crack. The curve for k^^(e) comes from the 
non $ i ngul ar part of the stress being removed from the crack 
by this iteration. The combination of the two designated 
k^(e) is what this iteration adds to the previous total 
stress intensity factor. (In this case, k(e) = k^(e) = 1.) 

Figure 3.9 
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Figure 3.10 


4. FUTURE APPLICATIONS - SEMI-ELLIPTICAL CRACK 

The problem of a flat elliptical crack lying completely 
within a half space and perpendicular to the surface has 
already been attempted in [9]. The results in [9] for the 
case in which the major axis of the ellipse lies in the sur- 
face (Figure 4.1) are labeled as approximate in recognition 
of the greater accuracy required due to the intersection of 
the crack front with the surface. The troublesome singular- 
ities encountered in the alternating method which were dis- 
cussed in detail in Section 3.4 are present in the semi- 
elliptical surface crack problem. Such considerations are 
neglected in [9], but other difficulties remain. The major 
difficulty is that the solution for the elliptical crack in 
an infinite medium is extremely complicated. 

Figure 4.1 

For the penny-shaped crack problem, the pressure on the 
crack needs only to permit a cosine series expansion of the 
form of equations (3.23) and (3.30). And the inclusion of 
more terms in the series is easily handled because of the 
nature of the solution. However, the solution for the ellip- 
tical crack requires the pressure to be given as a polynomial 
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(4.1) 


N M 

P(x.y) = l 1 a mn x m y n 

n=0 m=0 mn 

It is possible to work out the solution for any finite number 
of terms of equation (4.1), although, for a large number, it 
is impractical. Thus far, no way has been found to include 
the effect of additional terms in (4.1) by a recurrence scheme. 
The additional terms require the solution for the lower order 
terms and additional contributions combined in a complicated 
way. 


As a result of the difficulty, the elliptical crack prob- 
lem was restricted in [9] to be of the form 

p(x,y) = a Q0 + a 0 ,y + a 2Q x 2 + a Q2 y 2 + a 2 ,x 2 y + a Q3 y 3 

( 4 . 2 ) 

The six constants in this equation are determined by a least 
squares technique of fitting equation (4.2) to a number of 
known values of p(x,y). The odd powers of x are omitted 
because the problem is taken to be symmetric about the yz- 
plane. For the elliptical surface flaw, the pressure would 
also be made symmetric about the xz-plane, and the odd powers 
of y would be eliminated. Thus, only three constants, a Q0 , 
a 20* a 02’ wou ^ cl be ava il at) l e for the least square fitting 
necessary. Additional terms are clearly needed, but the 200 
terms used in the penny-shaped crack solution (40 in 0, 5 in r) 
cannot be approached. 


- 67 - 



As suggested by Kantorovich and Krylov [2], the Schwarz- 
Neumann Alternating Technique has wide applicability. It can 
be applied to other three-dimensional crack problems of 
interest such as the finite thickness crack problem. The 
crack may be taken to intersect the plate at any angle. Such 
problems would involve combined modes of crack extension even 
though the plate is loaded in simple tension. 

There is some arbitrariness in the choice of the two 
sequences of problems involved in the alternating method. For 
example, in Section 2.3, the stress, a(x), on the edge crack 
was extended symmetrically to negative values of x in the 
zeroth iteration. It would have been possible to make the 
stress skew-symmetric about the y-axis instead. The two se- 
quences of problems would then be quite different. Similarly, 
in the surface crack problem, the initial pressure on the 
penny-shaped crack could be taken to be skew-symmetric about 
the yz-plane. The problems which make up the two sequences 
would then have different characteristics. In general, the 
various alternatives should be studied to determine whether 
one provides a better representation of the essential proper- 
ties of the problem. 

It is conceivable that in some applications the alter- 
nating method will involve three or more sequences of solu- 
tions. In general, one sequence would be expected to be 
necessary for each part of the boundary of the material. 
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If there were three parts of the boundary and three solutions, 
the sequence would be formed similarly to the examples pre- 
sented here. Each solution would produce some undesirable 
residual stresses on the other two boundaries. Thus, each 
iteration would use each solution in turn to remove the 
residual stresses produced by the other two solutions in 
their previous applications. 
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APPENDIX 


Stress intensity factors for various particular loads 
on the edge crack are listed here. For other loads the 
general expression of eq . (2.50) may be used. In addition, 
results for the penny-shaped crack in an infinite body are 
given, includigg some recently published stress intensity 
factors for non-axi symmetri cal 1 y loaded cracks. 

Edge Cracks 

(1) Concentrated forces (Figure 2.4) 

For concentrated normal and shearing forces applied 
to the crack as shown in Figure 2.4, the stress intensity 
factors are 


= 2 Hfjb/a) 

71 /a*"- b' 2 


k 


2 



1+f (b/a) 
/a 2 - b 2 


where f(b/a) is given by equation (2.47) and Figure 2.5. 


Figure A-l 


(2) Linear stress (Figure A-l) 

In this case, the stress on the crack varies linearly 
from zero at the edge to twice the average stress at the 
crack tip (Figure A-l). The stress intensity factors in 
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terms of the average stresses are 


k = 1.3660 o /a 

1 a 

k, = 1 .3660 r /a 

Figure A- 2 

(3) Uniform stress near crack tip (Figure A-2) 

As shown in Figure A-2, the crack is free of stress 
on the part of its faces adjacent to the edge and is sub- 
jected to uniform stresses on the other part up to the crack 
tip. The stress intensity factors may be expressed as 

kj = £ arccos (b/a Jog/a [1 +h j ( b/a ) ] 

k 2 = | arccos(b/a)T 0 /a[l+h l (b/a)] 

where hj(b/a) is shown in Figure A-3. These results may 
also be expressed in terms of the average stresses on the 
crack , 

o g = ( 1 - b/ a ) a q x a = ( 1 -b/a 

The stress intensity factors 
kj = h 2 (b/a)o a /a 
k 2 = h 2 (b/a)x a /a 
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then show the expected increase when the same total load is 
applied nearer the crack tip. See Figure A-3, 

Figure A-3 

(4) Uniform stress (See Table I) 

In the case of uniform shear loading, the stress 
where the crack intersects the edge must be zero. See also 
the discussion in Section 2.4 of a similar restriction on 
the normal stress when the alternating method is used di- 
rectly. These restrictions do not affect the use of case 
(3) above for very small values of b/a. As b/a+0 it is found 
that 


kj = 1.1215 o 0 /a 
k, = 1 .1215 x n /a 

By superposition the first result may be applied to the case 
of uniform tension at infinity parallel to the free edge 
when the crack is stress free. There is no corresponding 
application of -the second result. 

(5) Polynomial stress [18] 

If the form of the stress in Figure A-l is not linear, 
but is given by 

o v (x,0) = -a Q l C (x/a) 
y u n = 0 
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where C 0 , •••, C 10 are arbitrary constants, then the 

stress intensity factor is [18] 

k = o 0 /7a[0.7930C 0 + 0.4829C 1 + 0 . 37 1 6C 2 + 0.3118C 3 

+ 0. 2 7 3 5 C 4 + 0 . 2464C 5 + 0.2260C g + 0.2099C, 

+ 0 . 1 9 6 8 C 0 + 0.1 858C 9 + 0.1 765C lfl ] 

If there is shear loading t ' (x,0) given by the expression 

xy 

for a (x,0), but with Oq replaced by Tq, then by the state- 
ment following eq.(2.50), the Mode II stress intensity 
factor, k , is given by the expression for k with. a n replaced 
by t q . 

Penny - Shaped Cracks 

(6) General result (Figure 3.1) 

The pressure applied to the crack is in the form 

oo 

a = -p 0 l K (r)cosme, 2=0 , r<l 
z m=0 

where r is the distance from the z-axis divided by the crack 
radius, a. The resulting Mode I stress intensity factor is 

ki = % P 0 /z l 9 m (l )cosme 
m=0 


where 


«.«*> 


0 


/I ~?2 
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These formulas apply only when the stress on the crack 
is symmetric about a plane containing the z-axis (the 8=0 
plane). For the special case of axisymmetric loading the 
only nonzero K m (r) is K Q (r), and the equations simplify 
consi derably . 

(7) Uniform tension (Figure A-4) [26] 

By superposition, the stress intensity factor due to 
the uniaxial tension may be determined by considering the 
same geometry, but with uniform pressure, 

°z = ~°0 ’ z=0, r<1 

on the crack and zero stress at infinity. For this problem 
the parameters of case (6) are 

p 0 = a 0 ’ K 0 (r) = 1 5 K m( r ) = °» m - 1 

Therefore, the stress intensity factor is [26]. 


Figure A-4 

(8) Uniform shear (Figure A-4) [30] 

The solution in [30] for constant shearing stresses 
on the faces of a penny-shaped crack may be used to obtain 
the stress intensity factors for Figure A-4. 
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1 4 y— 

k ! " ? T (/» sine 

k 3 = J T 0 /5 ' COse 

The uniform shear load produces no Mode I stress intensity 
factor. 


Figure A-5 


(9) Concentrated forces (Figure A-5) [6] 

The stress intensity factors for Modes II and III 
are zero, and 

oo 

k, = ~ /a — - — [^ + j r m cosm9 , cosm0] 

1 7T 2 n C. L «, 1 1 

ira 2 /I - r 2 m=l 

i 

The two pairs of concentrated forces act a distance of ar t 
out on the rays 0 = 10^ 


In the case of a single pair of forces of magnitude F on 
the crack, may be put equal to zero and 2P replaced by F 
in the above expression. If, further, the single pair of 
forces acts at the center of the crack. 
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(10) Linear stress [6] 

Case (6) may be applied to the case of normal stress 
varying linearly with x=rcos0. For 

o zz = -pprcosS, z = 0, r<l 

the stress varies from zero on the y-axis to a maximum of 
Pq at the points where the crack front intersects the x-axis 
of Figure 3.1. The only non-zero K m (r) is 

Kj ( r ) = r 

and the stress intensity factor, 

k i = I? po^ cos9 

is easily obtained. 

Figure A-6 

(11) Bending stress (Figure A-6 ) [63 

Cases (7) and (10) may be superposed to obtain the 
solution of the bending of a beam containing a small penny- 
shaped crack in a transverse section. If the crack is far 
enough from the lateral surfaces of the beam their inter- 
action may be neglected. Reference [6] specifies 

d > 3a c-b > 3a 
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where the dimensions of the beam are shown in Figure A-6. 

An additional restriction must be placed on the dimensions 
to insure that the crack faces are not compressed together, 
for the solution would be invalid if that happened. This 
restriction, 

b > I a 

simply makes the stress intensity factor positive at all 
points of the crack front. In terms of the moment of inertia, 
I = 2dc 3 /3, of the cross-section, the stress produced in an 
uncracked beam is 

M Q a 

o zz = — — (b/a - rcose), z = 0 , r<l 

By superposition, the stress intensity factor 
? M n a 

ki = - — j— /a (b/a - (2/3)cos6) 

may be found by considering the beam to be loaded only by 
stresses on the faces of the crack equal and opposite to those 
produced in the uncracked beam. 

Surface Crack 

(12) Uniform stress 

The geometry of the half penny-shaped crack at the 
surface of a half space is shown in Figure 3.6. For the case 
of uniform tension at infinity perpendicular to the crack 
(or uniform pressure on the crack faces), the stress inten- 
sity factor is given in Figure 3.8. 
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Fig. 2.1: Edge-cracked plate 
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Fig. 2. 



2: Cracked infinite plate 
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Fig. 2.3: Edge-loaded half plane 
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Fig. 3.4: Auxiliary coordinates for singular solutions 
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Fig. A - 1 : Linearly loaded edge crack 
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Fig. A - 2 : Partially loaded edge crack 
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